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Abstract 

A non-technical description of the Operator Product Expansion and Renormalization Group 
techniques as applied to weak decays of mesons is presented. We use this opportunity to 
summarize briefly the present status of the next-to-leading QCD corrections to weak decays 
and their implications for the unitarity triangle, the ratio e' /e, the radiative decay B — > Xg'j, 
and the rare decays K + — > Tt + vv and Kl — > -k^vv. 
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1 Preface 



It is a great privilege and a great pleasure to give this talk at the symposium celebrating 
the 70th birthday of Wolfhart Zimmermann. The Operator Product Expansion Q to which 
Wolfhart Zimmermann contributed in such an important manner B ||, |[ had an important 
impact on my research during the last 20 years. I do hope very much to give another talk 
on this subject in 2008 at a symposium celebrating Wolfhart Zimmermanns 80th birthday. I 
am convinced that OPE will play an important role in the next 10 years in the field of weak 
decays as it played already in almost 25 years since the pioneering applications of this very 
powerful method by Gaillard and Lee ||] and Altarelli and Maiani ||. 



2 Operator Product Expansion 

The basic starting point for any serious phenomenology of weak decays of hadrons is the 
effective weak Hamiltonian which has the following generic structure 

^// = ^E^ckmCU/^ • (2-1) 

Here Gf is the Fermi constant and Qi are the relevant local operators which govern the 
decays in question. The Cabibbo-Kobayashi-Maskawa factors Vq KM |7|, || and the Wilson 
Coefficients C, (lj describe the strength with which a given operator enters the Hamiltonian. 
In the simplest case of the /3-decay, 7~C e ff takes the familiar form 

n[f f = -| cos0 c [« 7m (1 - 75 )d ® erf {I - 75 K] , (2.2) 

where V u d has been expressed in terms of the Cabibbo angle. In this particular case the Wilson 
Coefficient is equal unity and the local operator, the object between the square brackets, is 
given by a product of two V — A currents. This local operator is represented by the diagram 
(b) in fig. |l|. Equation (^^) represents the Fermi theory for /3-decays as formulated by 
Sudarshan and Marshak [|| and Feynman and Gell-Mann |]l0| forty years ago, except that 
in (^) the quark language has been used and following Cabibbo a small departure of V u d 
from unity has been incorporated. In this context the basic formula (^^) can be regarded 
as a generalization of the Fermi Theory to include all known quarks and leptons as well as 
their strong and electroweak interactions as summarized by the Standard Model. It should 
be stressed that the formulation of weak decays in terms of effective Hamiltonians is very 
suitable for the inclusion of new physics effects. We will discuss this issue briefly later on. 
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(a) (b) 



Figure 1: /3-decay at the quark level in the full (a) and effective (b) theory. 

Now, I am aware of the fact that the formal operator language used here is hated by 
experimentalists and frequently disliked by more phenomenological minded theorists. Con- 



sequently the literature on weak decays, in particular on B-meson decays [11|, is governed 
by Feynman diagram drawings with W-, Z- and top quark exchanges, rather than by the 
operators in In the case of the /3-decay we have the diagram (a) in fig. |l[ Yet such 

Feynman diagrams with full W-propagators, Z-propagators and top-quark propagators really 
represent the situation at very short distance scales 0(M\y,z> rnt), whereas the true picture 
of a decaying hadron with masses 0{m^, m c , rax) is more properly described by effective 
point-like vertices which are represented by the local operators Qi. The Wilson coefficients 
Ci can then be regarded as coupling constants associated with these effective vertices. 



Thus 7~ieff m (2-1) is simply a series of effective vertices multiplied by effective coupling 



constants Cj. This series is known under the name of the operator product expansion (OPE) 



[H]~0' [12]- Due to the interplay of electroweak and strong interactions the structure of 
the local operators (vertices) is much richer than in the case of the /3-decay. They can be 
classified with respect to the Dirac structure, colour structure and the type of quarks and 
leptons relevant for a given decay. Of particular interest are the operators involving quarks 
only. They govern the non-leptonic decays. To be specific let us list the operators relevant 
for non-leptonic B-meson decays. They are: 
Current Current : 

Ql = (c a bj3) V -A (S[3C a ) V -A Q2 = (cb)v~A (sc)v-A (2-3) 
QCD Penguins : 

Q3 = {sb) V -A (Q9)v-A Qi = {sJ>p)v-A (qpQa)v-A (2.4) 

q=u,d,s,c,b q=u,d,s,c,b 
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Q 5 = {sb)v- A Y (.QlW+A Qg = {s a bf3)v-A Y (Q0a)v+A 

q=u,d,s,c,b 



q=u,d,s,c,b 

Electroweak Penguins : 

3 

Qt = t. (sb) V -A Y e i (QQ)v+A 



Qs 



q=u,d,s,c,b 



{s a bf3) V -A Y e qi<lPla)v+A 
q=u,d,s,c,b 



Qg 



- {sb) V -A Y e <?(<?<?)v-A QlQ = ~ (s a bf3)v-A Y e 1 (WQodv-A 
~ q=u,d,s,c,b ~ q=u,d,s,c,b 



(2.5) 



(2.6) 



(2.7) 



Here, a and (5 are colour indices and e q denotes the electrical quark charges reflecting the 
electroweak origin of Q-j, . . . , (Jio- Q2, Q3-6 and Qj : g originate in the tree level iy ± -exchange, 
gluon penguin and (7, Z°)-penguin diagrams respectively. These are the diagrams a)— c) in 
fig. |2[ To generate Q\, Qg and Q\q additional gluonic exchanges are needed. The operators 
given above have dimension six. Of interest are also operators of dimension five which are 
responsible for the B — > S7 decay. They originate in the diagram d) in fig. || where 7 and 
the gluon are on-shell. They will be given in Section 7. In what follows we will neglect the 
higher dimensional operators as their contributions to weak decays are marginal. 



u,c,t 



u,c,t 




Figure 2: Typical Penguin and Box Diagrams. 
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Now what about the couplings Cj(/i) and the scale /U? The important point is that 
Ci(fi) summarize the physics contributions from scales higher than [i and due to asymptotic 
freedom of QCD they can be calculated in perturbation theory as long as fi is not too small. 
Ci include the top quark contributions and contributions from other heavy particles such as 
W, Z-bosons and charged Higgs particles or supersymmetric particles in the supersymmetric 
extensions of the Standard Model. At higher orders in the electroweak coupling the neutral 
Higgs may also contribute. Consequently Cj(yu) depend generally on mt and also on the 
masses of new particles if extensions of the Standard Model are considered. This dependence 
can be found by evaluating the box and penguin diagrams with full W-, Z-, top- and new 
particles exchanges shown in fig. |2| and properly including short distance QCD effects. The 
latter govern the ^-dependence of the couplings Ci(/i). 

The value of fi can be chosen arbitrarily. It serves to separate the physics contributions 
to a given decay amplitude into short-distance contributions at scales higher than \x and 
long-distance contributions corresponding to scales lower than [i. It is customary to choose [i 
to be of the order of the mass of the decaying hadron. This is 0{m\ ) ) and 0{m c ) for B-decays 
and D-decays respectively. In the case of K-decays the typical choice is \x = 0(1 — 2 GeV) 
instead of 0{jxik\ which is much too low for any perturbative calculation of the couplings 
Ci. 

Now due to the fact that \i <C Mw,Zi m t, large logarithms lnA%/^ compensate in the 
evaluation of Cj(/z) the smallness of the QCD coupling constant a s and terms a™ (In M\y/V) n , 
a" (In M\v / ^) n ~ 1 etc. have to be resummed to all orders in a s before a reliable result for Cj can 
be obtained. This can be done very efficiently by means of the renormalization group methods 
fll3| , |l4| , 15]. Indeed solving the renormalization group equations for the Wilson coefficients 
Cj(/i) summs automatically large logarithms. The resulting renormalization group improved 
perturbative expansion for Cj(/x) in terms of the effective coupling constant a s (n) does not 
involve large logarithms and is more reliable. 

It should be stressed at this point that the construction of the effective Hamiltonian H. e ff 
by means of the operator product expansion and the renormalization group methods can be 
done fully in the perturbative framework. The fact that the decaying hadrons are bound 
states of quarks is irrelevant for this construction. Consequently the coefficients Ci(ji) are 
independent of the particular decay considered in the same manner in which the usual gauge 
couplings are universal and process independent. 

Having constructed the effective Hamiltonian we can proceed to evaluate the decay 
amplitudes. An amplitude for a decay of a given meson M = K,B,.. into a final state 
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F = ixvv, 7T7T, DK is simply given by 

A (M - F) = (F\n eff \M) = ^J2 V CKMCi^)(F\Q i ^)\M), (2.8) 

where (F\Qi(jj)\M) are the hadronic matrix elements of Qi between M and F. As indicated 
in fl2.8|) these matrix elements depend similarly to Cj(/i) on [i. They summarize the physics 
contributions to the amplitude A(M — * F) from scales lower than fj,. 

We realize now the essential virtue of OPE: it allows to separate the problem of calcu- 
lating the amplitude A(M — > F) into two distinct parts: the short distance (perturbative) 
calculation of the couplings Cj(/z) and the long-distance (generally non-perturbative) calcu- 
lation of the matrix elements (Qi (//)). The scale /i, as advertised above, separates then the 
physics contributions into short distance contributions contained in Cj(/x) and the long dis- 
tance contributions contained in (Qi(/j,)). By evolving this scale from \x = 0(M\y) down to 
lower values one simply transforms the physics contributions at scales higher than /i from the 
hadronic matrix elements into Cj(/i). Since no information is lost this way the full amplitude 
cannot depend on \x. Therefore the /i-dependence of the couplings Cj(/x) has to cancel the //- 
dependence of (QiGu)). In other words it is a matter of choice what exactly belongs to Ci(fi) 
and what to (Qj (//))• This cancellation of /i-dependence involves generally several terms in 
the expansion in ( |2,8| ). 

Clearly, in order to calculate the amplitude A(M — > F), the matrix elements (Qi(/i)) have 
to be evaluated. Since they involve long distance contributions one is forced in this case to use 
non-perturbative methods such as lattice calculations, the 1/N expansion (N is the number 
of colours), QCD sum rules, hadronic sum rules, chiral perturbation theory and so on. In 
the case of certain B-meson decays, the Heavy Quark Effective Theory (HQET) turns out to 
be a useful tool. Needless to say, all these non-perturbative methods have some limitations. 
Consequently the dominant theoretical uncertainties in the decay amplitudes reside in the 
matrix elements (Qj (//)). 

The fact that in most cases the matrix elements (Qi(fJ-)) cannot be reliably calculated at 
present, is very unfortunate. One of the main goals of the experimental studies of weak decays 
is the determination of the CKM factors V^ KM and the search for the physics beyond the 
Standard Model. Without a reliable estimate of (Qi(/i)) this goal cannot be achieved unless 
these matrix elements can be determined experimentally or removed from the final measurable 
quantities by taking the ratios or suitable combinations of amplitudes or branching ratios. 
However, this can be achieved only in a handful of decays and generally one has to face 
directly the calculation of (Qi(fi)). 

Now in the case of semi-leptonic decays, in which there is at most one hadron in the 
final state, the chiral perturbation theory in the case of K-decays and HQET in the case of 
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B-decays have already provided useful estimates of the relevant matrix elements. This way 
it was possible to achieve satisfactory determinations of the CKM elements V us and in 
K — * Tiev and B — > D*ev respectively. Similarly certain rare decays like K — > nvv and 
B — > /i/2 can be calculated very reliably. 

The case of non-leptonic decays in which the final state consists exclusively out of hadrons 
is a completely different story. Here even the matrix elements entering the simplest decays, 
the two-body decays like K — > tttt, D — > Kir or B — > DK cannot be calculated in QCD 
reliably at present. More promising in this respect is the evaluation of hadronic matrix 
elements relevant for K° — K° and B® 3 — B® s mixings. 

Returning to the Wilson coefficients Cj(/u) it should be stressed that similar to the effective 
coupling constants they do not depend only on the scale \x but also on the renormalization 
scheme used: this time on the scheme for the renormalization of local operators. That the 
local operators undergo renormalization is not surprising. After all they represent effective 
vertices and as the usual vertices in a field theory they have to be renormalized when quan- 
tum corrections like QCD or QED corrections are taken into account. As a consequence of 
this, the hadronic matrix elements {Qi(fi)) are renormalization scheme dependent and this 
scheme dependence must be cancelled by the one of Cj(^) so that the physical amplitudes 
are renormalization scheme independent. Again, as in the case of the ^-dependence, the can- 
cellation of the renormalization scheme dependence involves generally several terms in the 



expansion (2.8). 

Now the [i and the renormalization scheme dependences of the couplings Cj(/i) can be eval- 
uated efficiently in the renormalization group improved perturbation theory. Unfortunately 
the incorporation of these dependences in the non-perturbative evaluation of the matrix ele- 
ments {Qi(n)) remains as an important challenge and most of the non-perturbative methods 
on the market are insensitive to these dependences. The consequence of this unfortunate 
situation is obvious: the resulting decay amplitudes are \x and renormalization scheme de- 
pendent which introduces potential theoretical uncertainty in the predictions. On the other 
hand in certain decays these dependences can be put under control. 

So far I have discussed only exclusive decays. It turns out that in the case of inclusive 
decays of heavy mesons, like B-mesons, things turn out to be easier. In an inclusive decay 
one sums over all (or over a special class) of accessible final states so that the amplitude for 
an inclusive decay takes the form: 

A(B - X) = ^ J2 V^ KM C^){f\Q^)\B) . (2.9) 

V 1 feX 

At first sight things look as complicated as in the case of exclusive decays. It turns out, 
however, that the resulting branching ratio can be calculated in the expansion in inverse 
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powers of mb with the leading term described by the spectator model in which the B-meson 
decay is modelled by the decay of the 6-quark: 

Br(S -► X) = Br(6 -> q) + C7(— y) . (2.10) 

m h 

This formula is known under the name of the Heavy Quark Expansion (HQE) J16||-[18|. Since 
the leading term in this expansion represents the decay of the quark, it can be calculated in 
perturbation theory or more correctly in the renormalization group improved perturbation 
theory. It should be realized that also here the basic starting point is the effective Hamiltonian 



(2.1) and that the knowledge of the couplings Cj(/i) is essential for the evaluation of the 
leading term in ( |2.10| ). But there is an important difference relative to the exclusive case: the 
matrix elements of the operators Qi can be "effectively" evaluated in perturbation theory. 
This means, in particular, that their fi and renormalization scheme dependences can be 
evaluated and the cancellation of these dependences by those present in Ci(p) can be explicitly 
investigated. 

Clearly in order to complete the evaluation of Br(B — > X) also the remaining terms in 
( 2.10| ) have to be considered. These terms are of a non-perturbative origin, but fortunately 
they are suppressed by at least two powers of m&. They have been studied by several authors 
in the literature with the result that they affect various branching ratios by less than 10% 
and often by only a few percent. Consequently the inclusive decays give generally more 
precise theoretical predictions at present than the exclusive decays. On the other hand their 
measurements are harder. There is of course an important theoretical issue related to the 
validity of HQE in (plo| ) which appear in the literature under the name of quark-hadron 
duality. I will not discuss it here. Recent discussions of this issue can be found in [jl9[ . 

We have learned now that the matrix elements of Qi are easier to handle in inclusive decays 
than in the exclusive ones. On the other hand the evaluation of the couplings Ci(fi) is equally 
difficult in both cases although as stated above it can be done in a perturbative framework. 
Still in order to achieve sufficient precision for the theoretical predictions it is desirable to 
have accurate values of these couplings. Indeed it has been realized at the end of the eighties 
that the leading term (LO) in the renormalization group improved perturbation theory, in 
which the terms a"(lnM\v//j) n are summed, is generally insufficient and the inclusion of 
next-to- leading corrections (NLO) which correspond to summing the terms a"(lnil%//i) Il ~ 1 
is necessary. In particular, unphysical left-over /^-dependences in the decay amplitudes and 
branching ratios resulting from the truncation of the perturbative series are considerably 
reduced by including NLO corrections. These corrections are known by now for the most 
important and interesting decays and will be briefly reviewed below. 
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3 Penguin— Box Expansion and OPE 



The FCNC decays, in particular rare and CP violating decays are governed by various penguin 
and box diagrams with internal top quark and charm quark exchanges. Some examples are 
shown in fig. 2. These diagrams can be evaluated in the full theory and are summarized 
by a set of basic universal (process independent) m t -dependent functions F r (xt) pO] where 



xt = mi /M w . Explicit expressions for these functions can be found in |2l| , p?2j , |23[| . 

It is useful to express the OPE formula ( |2.g| ) directly in terms of the functions F r (xt) (2 
To this end we rewrite the A{M — > F) in (|2.8|) as follows 



A (M - F) = %^ckm J2( F I °M I M ) Uki (M. M w ) Q(M w ), (3.11) 

where Ukj(fJ>, M\y) is the renormalization group transformation from Mw down to /U. Explicit 
formula for this transformation will be given below. In order to simplify the presentation we 
have removed the index "i" from V^. KM 

Now Cj(M\v) are linear combinations of the basic functions F r (xt) so that we can write 

d(M w ) = Ci + Y, h ir F r (x t ) (3.12) 



where Cj and hi r are m t -independent constants. Inserting ( p.!2[ ) into (3.11) and summing 
over % and k we find 

A (M — > F) = P (M ~^F)+J2 Pr{M F) F r (x t ), (3.13) 

r 

with 



P (M — > F) = ^(F | O k (fi) | M) C> fci (/x, M W )Q , (3.14) 
P r ( M — > F) = ^ (F | O fc (/i) | M) U kl {fi,M w )h ir , (3.15) 

i,k 

where we have suppressed the overall factor (Gp /y2)VcKM- I would like to call fl3,13| ) 
Penguin-Box Expansion (PBE) p5| . 

The coefficients Pq and P r are process dependent. This process dependence enters through 
(F | Ok(n) | M). In certain cases like K — > ttvu these matrix elements are very simple 
implying simple formulae for the coefficients Pq and P r . In other situations, like e' /e, this is 
not the case. 

Originally PBE was designed to expose the mt-dependence of FCNC processes p5|| . After 
the top quark mass has been measured precisely this role of PBE is less important. On the 
other hand, PBE is very well suited for the study of the extentions of the Standard Model 
in which new particles are exchanged in the loops. We know already that these particles are 
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heavier than W-bosons and consequently they can be integrated out together with the weak 
bosons and the top quark. If there are no new local operators the mere change is to modify 
the functions F r {xt) which now acquire the dependence on the masses of new particles such 
as charged Higgs particles and supersymmetric particles. The process dependent coefficients 
Pq and P r remain unchanged. This is particularly useful as the most difficult part is the 
evaluation of Ukj(n, Myy) and of the hadronic matrix elements, both contained in these 
coefficients. However, if new effective operators with different Dirac and colour structures 
are present the values of Po an d P r are modified. Examples of the applications of PBE to 
physics beyond the Standard Model can be found in |2(| ^7], [2gfl . 

The universality of the functions F r (xt) can be violated partly when QCD corrections to 
one loop penguin and box diagrams are included. For instance in the case of semi-leptonic 
FCNC transitions there is no gluon exchange in a Z°-penguin diagram parallel to the Z°- 
propagator but such an exchange takes place in non-leptonic decays in which the bottom line 
is a quark-line. Thus the general universality of F r {xt) present at one loop level is reduced 
to two universality classes relevant for semi-leptonic and non-leptonic transitions. However, 
the 0(a s ) corrections to the functions F r {xt) are generally rather small when the top quark 
mass rn t (m t ) is used and consequently the inclusion of QCD effects plays mainly the role in 
reducing various /i-dependences. 

In order to see the general structure of A(M —* F) more transparently let us write it as 
follows: 

A(M — F) = B M -+FV C KMVQCDF(x t ) + Charm (3.16) 

where the first term represents the internal top quark contribution and "Charm" stands for 
remaining contributions, in particular those with internal charm quark exchanges. F{xt) 
represents one of the universal functions and ??qcd the corresponding short distance QCD 
corrections. The parameter Bm^f represents the relevant hadronic matrix element, which 
can only be calculated by means of non-perturbative methods. However, in certain lucky 
situations Bm-*f can be extracted from well measured leading decays and when it enters 
also other decays, the latter are then free from hadronic uncertainties and offer very useful 
means for extraction of CKM parameters. One such example is the decay K + — * ir + vi) for 
which one has 



Br{K + vT 



vv) 



o 



QED 



Br{K" 



K? s 2tt 2 S i n 4 e w 



V^Vu^Fixt) + V: s V cdV c QCD F(a 



(3.17) 

The factor in square brackets stands for the "B — factor" in ( |3.16 ), which is given in terms 
of well measured quantities. Since V cs , V c d and Vt s are already rather well determined and 
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F(xi) and f/o CD can be calculated in perturbation theory, the element Vtd can be extracted 
from Br(K + — > tt + ui>) without essentially any theoretical uncertainties. We will be more 
specific about this in Section 7. 

4 Motivations for NLO Calculations 

Going beyond the LO approximation for Ci(fi) is certainly an important but a non-trivial 
step. For this reason one needs some motivations to perform this step. Here are the main 
reasons for going beyond LO: 

• The NLO is first of all necessary to test the validity of the renormalization group 
improved perturbation theory. 

• Without going to NLO the QCD scale A-j^g |2^] extracted from various high energy 
processes cannot be used meaningfully in weak decays. 

• Due to renormalization group invariance the physical amplitudes do not depend on the 
scales fi present in a s or in the running quark masses, in particular mt(fx), mb(/i) and 
m c ((j). However, in perturbation theory this property is broken through the truncation 
of the perturbative series. Consequently one finds sizable scale ambiguities in the 
leading order, which can be reduced considerably by going to NLO. 

• The Wilson Coefficients are renormalization scheme dependent quantities. This scheme 
dependence appears first at NLO. For a proper matching of the short distance con- 
tributions to the long distance matrix elements obtained from lattice calculations it is 
essential to calculate NLO. The same is true for inclusive heavy quark decays in which 
the hadron decay can be modeled by a decay of a heavy quark and the matrix elements 
of Qi can be effectively calculated in an expansion in l/m^. 

• In several cases the central issue of the top quark mass dependence is strictly a NLO 
effect. 

5 General Structure of Wilson Coefficients 

We will give here a formula for the Wilson coefficient C(/i) of a single operator Q including 
NLO corrections. The case of several operators which mix under renormalization is much 
more complicated. Explicit formulae are given in |2l] , |23| . 
C(p) is given by 

C(p) = U(ji, M W )C(M W ) (5.18) 
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where 



U(fj,, Mw) = exp 



9S(P) d g nQ{9 ' s) 

g s (M w ) P(g's) 



(5.19) 



is the evolution function, which allows to calculate C(/i) once C(Myj) is known. The latter 
can be calculated in perturbation theory in the process of integrating out W^, Z° and top 



quark fields. Details can be found in [21, 23]. Next is the anomalous dimension of the 
operator Q and (3{g s ) is the renormalization group function which governs the evolution of 
the QCD coupling constant 
At NLO we have 

ct s {M w ) 



1 + 



4tt 



-B 



7Q(« S ) =1q T-+7q I — I 



-Po 



47T 

16vr 2 



ft 



4-7T / 



(16vr 2 ) 2 



Inserting the last two formulae into ( 5.19] ) and expanding in a s we find 



U(n,M w ) 



i + ^Mj 



47T 



a s (M ] 



w 



a s {M w ) 

47T 



J 



with 



J = J A 



-7 (1) 



(o) 
7q 



(5.20) 
(5.21) 
(5.22) 

(5.23) 
(5.24) 



2/? 2/3 

Inserting ( ggD and (^20]) into ( {Up we find an important formula for C(/i) in the NLO 
approximation : 

d 



C(/i) 



1 + ^-^J 

47T 



a s (M 



w 



1 + ^^1(B-J) 
4ir 



(5.25) 



6 Status of NLO Calculations 

Since the pioneering leading order calculations of Wilson coefficients for current -current |5], ^| 
and penguin operators pC[], enormous progress has been made, so that at present most of the 
decay amplitudes are known at the NLO level. We list all existing NLO calculations for weak 
decays in table |l[ In addition to the calculations in the Standard Model we list the calculations 
in two-Higgs doublet models and super symmetry. In table ^ we list references to calculations 
of two- loop electroweak contributions to rare decays. The latter calculations allow to reduce 
scheme and scale dependences related to the definition of electroweak parameters like sin 2 6yy , 



cxqed, etc. Next, useful techniques for three-loop calculations can be found in [77] and a very 
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general discussion of the evanescent operators including earlier references is presented in 78]. 
Further details on these calculations can be found in the orignal papers, in the review [21] 
and in the Les Houches lectures [ p3[| . Some of the implications of these calculations will be 
analyzed briefly in subsequent sections. 

7 Applications: News 
7.1 Preliminaries 

There is a vast literature on the applications of NLO calculations listed in table 1. As they 
are already reviewed in detail in [^l], [2^, [2*|] there is no point to review them here again. I 
will rather discuss briefly some of the most important applications in general terms. This 



will also give me the opportunity to update some of the numerical results presented in 23]. 
This update is related mainly to the improved experimental lower bound on B® — B® mixing 
((AM)j > 12.4/ps) and a slight increase in (V^I/jV^I: 0.091 ±0.016, both presented at the 
last Rochester Conference in Vancouver Iff' 

A=(p,f|) 




p-KTi / \. l-p4r| 



C=(0,0) B=(1,0) 

Figure 3: Unitarity Triangle. 



7.2 Unitarity Triangle 

The standard analysis of the unitarity triangle (see fig. ||) uses the values of | V^ s | , \V c b\, 
\V u b/Vcb\ extracted from tree level K- and B- decays, the indirect CP-violation in Ki — * tttt 
represented by the parameter e and the B® — B® mixings described by the mass differences 
(AM)ds- From this analysis follows the allowed range for (g,fj) describing the apex of the 
unitarity triangle. Here [^J] 

e=e(l-y), V = r,(l-¥-). (7.26) 
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where A, g and rj are Wolfenstein parameters [82] with \ V US \ = A 
V ub = X\V cb \(g-i7 1 ), V td = X\V cb \(l - 



0.22. We have in particular 
-if]). (7.27) 



T) 7^ is responsible for CP violation in the Standard Model. 

The allowed region for (g, ff) is presented in fig. It is the shaded area on the right 
hand side of the solid circle which represents the upper bound for (AM)^/(AM) S . The 
hyperbolas give the constraint from e and the two circles centered at (0, 0) the constraint 
from The white areas between the lower e- hyperbola and the shaded region are 

excluded by B® — B® mixing. We observe that the region g < is practically excluded. The 
main remaining theoretical uncertainties in this analysis are the values of non-perturbative 
parameters: B K in e, F Bd ^W d in (AM) d and £ = F Bs ^W s / F Bd ^W d in (AM) d /(AM) s . I 
have used B K = 0.80 ± 0.15, F Bd ^fB~ d = 200 ± 40MeV and £ max = 1.2. On the experimental 
side [V^b/Vcbl and (AM) S should be improved. 



0.4 



0.0 



K,(AM),).(1.2,12.4ps 




-0.4 -0.2 0.0 0.2 0.4 0.6 

P 



Figure 4: Unitarity Triangle 1998. 
From this analysis we extract 

\V U \ = (8.6 ± 1.6) • 10- 3 , lm(V t *V td ) = (1.38 ± 0.33) • 10~ 4 



and 



7.3 e'/e 



sin 2/5 = 0.71 ±0.13, sin 7 = 0.83 ± 0.17 



(7.28) 
(7.29) 



e'/e is the ratio of the direct and indirect CP violation in Kl — * vi"7r. A measurement of 
a non- vanishing value of e'/e would give the first signal for the direct CP violation ruling 
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out the superweak models [83]. In the Standard Model e'/e is governed by QCD penguins 
and electroweak (EW) penguins. The corresponding operators are given in ( |2 . 4[ ) - (pT7|) . With 
increasing m t the EW penguins become increasingly important 84, 85], and entering e'/e 
with the opposite sign to QCD penguins suppress this ratio for large m t . For m t ps 200 GeV 
the ratio can even be zero [85] ] . Because of this strong cancellation between two dominant 
contributions and due to uncertainties related to hadronic matrix elements of the relevant 
local operators, a precise prediction of e'/e is not possible at present. 

A very simplified formula (not to be used for any serious numerical analysis) which exhibits 
main uncertainties is given as follows 

0.8 



- = 15 • 10- 4 

e 



1.3 • 10- 4 



120 MeV 



m s (2 GeV) 



MS 



300 MeV 



[B 6 - Z(x t )B 8 



(7.30) 



where Z(xt) ~ 0.18(m t /M\v) 1 ' 86 represents the leading m t -dependence of EW penguins. 
Bq and i?8 represent the hadronic matrix elements of the dominant QCD-penguin operator 

)) respectively. 



Qq and the dominant electroweak penguin operator Q$ (see (2J5) and (£ 
Together with m s (2GeV) the values of these parameters constitute the main theoretical 



uncertainty in evaluating e'/e. Present status of m s , Bq and B$ is reviewed in [86, 23]. 
Roughly one has B§ = 1.0 ± 0.2 and B$ = 0.7 ± 0.2. Taking these values, r\ of fig. || and 
(V^l = 0.040 ±0.003, I find: 



e'/e 



(7.31) 



(5.7 ± 3.6) • 10~ 4 , m s (2 GeV) = 130 ± 20 MeV 
(9.1 ± 5.7) • 10~ 4 , m s (2 GeV) = 110 ± 20 MeV. 

where the chosen values for m s are in the ball park of various QCD sum rules and lattice 
estimates j86|. This should be compared with the result of NA31 collaboration at CERN 
which finds (e'/e) = (23 ± 7) • 10~ 4 |7| and the value of E731 at Fermilab, (e'/e) = (7.4 ± 
5.9) • 10~ 4 @. 

The Standard Model expectations are closer to the Fermilab result, but due to large 
theoretical and experimental errors no firm conclusion can be reached at present. The new 
improved data from CERN and Fermilab in 1999 and later from DA^NE should shed more 
light on e'/e. In this context improved estimates of Bq, B% and m s are clearly desirable. 



7.4 B -> X s7 

A lot of efforts have been put into predicting the branching ratio for the inclusive decay 
B — ► Xs'y including NLO QCD corrections and higher order electroweak corrections. The 



relevant references are given in table 1 and in [23], where details can be found. The final 
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result of these efforts can be summarized by 

Br(B -> X s j) th = (3.30 ± 0.15(scale) ± 0.26(par)) • 10" 4 (7.32) 

where the first error represents residual scale dependences and the second error is due to 
uncertainties in input parameters. The main achievement is the reduction of the scale de- 



pendence through NLO calculations, in particular those given in [61] and In the leading 
order the corresponding error would be roughly ±0.6 |${], 9C]. 



The theoretical result in ( [7.32 ) should be compared with experimental dat 



a: 



(3.15 ± 0.35 ± 0.41) • 10" 4 , CLEO 
Br(B -» X sl ex P = { ) / 7 - 33 

[ (3.11 ±0.80 ±0.72) • 10~ 4 , ALEPH, 

which implies the combined branching ratio: 

Br(B X s7 ) cxp = (3.14 ± 0.48) • 10~ 4 (7.34) 

Clearly, the Standard Model result agrees well with the data. In order to see whether any 
new physics can be seen in this decay, the theoretical and in particular experimental errors 
should be reduced. This is certainly a very difficult task. Most recent analyses of B — > X s j 
in supersymmetric models and two-Higgs doublet models are listed in table ||. 

7.5 K L — > ix Q vu and K + — > ix + vv 

Kl — ► iftvv and K + — > -k + vv are the theoretically cleanest decays in the field of rare K- 
decays. Kl — ► ix^vv is dominated by short distance loop diagrams (Z-penguins and box 
diagrams) involving the top quark. K + — * ix^vv receives additional sizable contributions 
from internal charm exchanges. The great virtue of Kl — > -k^vv is that it proceeds almost 



exclusively through direct CP violation [91| and as such is the cleanest decay to measure this 
important phenomenon. It also offers a clean determination of the Wolfenstein parameter 
rj and in particular offers the cleanest measurement of ImV^ s Vtd [92]. K + — ► 7r + i/i> is CP 



conserving and offers a clean determination of \Vtd\- Due to the presence of the charm 
contribution and the related m c dependence it has a small scale uncertainty absent in Kl — > 



The next-to-leading QCD corrections |49|, |50j, |53|, |51], |52|] to both decays considerably 
reduced the theoretical uncertainty due to the choice of the renormalization scales present 
in the leading order expressions, in particular in the charm contribution to K + — > tt + vu. 
Since the relevant hadronic matrix elements of the weak currents entering K — > ttvu can be 
related using isospin symmetry to the leading decay K + — ► the resulting theoretical 

expressions for Br( Kl — ► n^vv) and Br(i^T + — > ^uv) are only functions of the CKM 
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parameters, the QCD scale A^jg and the quark masses mt and m c . The isospin braking 



corrections have been calculated in [93]. The long distance contributions to K + — > n + vv have 
been considered in |)4j and found to be very small: a few percent of the charm contribution 
to the amplitude at most, which is safely neglegible. The long distance contributions to 
Kl — ► iflvv are negligible as well p5| . 

The explicit expressions for Br(K + — > 7r + i/p) and Br(Ki — > -k^uv) can be found in 



1 21, E2, Here we give approximate expressions in order to exhibit various dependences: 



Br(K + -> 7r + i/z?) = 0.7 • 10 



-10 



w. 



id 



0.010 



r6 



0.040 



m t {m t ) 



170 GeV 



2.3 



+ cc + tc 



Br{K L 



3.0 • 10 



-ii 



0.39 



m t (m t ) 



170 GeV 



2.3 



V, 



cb 



0.040 



(7.35) 



(7.36) 



where in (7.35) we have shown explicitly only the pure top contribution. 

The impact of NLO calculations is the reduction of scale uncertainties in Br{K + — > n + vV) 
from ±23% to ±7%. This corresponds to the reduction in the uncertainty in the determination 
of \Vtd\ from ±14% to ±4%. The remaining scale uncertainties in Br(Ki — > tt^uv) and in 
the determination of r\ are fully negligible. 

Updating the analysis of [03] one finds 52]: 



Br(K A 



it'vv) 



1.2 ±3.2) • 10- 



-ii 



Br(K L -» n Q vv) = (3.1 ± 1.3) • 10 



-ii 



(7.37) 



where the errors come dominantly from the uncertainties in the CKM parameters. 



As stressed in [92] simultaneous measurements of Br(K + — ► tt + vv) and Br{K^ — > ir v&) 
should allow a clean determination of the unitarity triangle as shown in fig. ||. In particu- 
lar the measurements of these branching ratios with an error of ±10% will determine \Vtd\, 
lmV t *V t d and sin 2/3 with an accuraccy of ±10%, ±5% and ±0.05 respectively. The compari- 
sion of this determination of sin 2/3 with the one by means of the CP-asymmetry in B — > ipKs 
should offer a very good test of the Standard Model. 

Experimentally we have [p6]| 



BR(K + -> tt + vu) = (4.2+9.7-3.5) • 10 



-10 



(7.38) 



and the bound 



BR{K L -» t^vv) < 1.6 • 10~ b . (7.39) 
Moreover from ( |7.38| ) and isospin symmetry one has |)8| BR(Kl — ► ifivv) < 6.1 • 10 -9 . 



The central value in (7.38) is by a factor of 4 above the Standard Model expectation 
but in view of large errors the result is compatible with the Standard Model. The analysis 
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(0,0) 



(1,0) 



Figure 5: Unitarity triangle from K 



■kvv. 



of additional data on K + — * ix^vv present on tape at BNL787 should narrow this range in 
the near future considerably. In view of the clean character of this decay a measurement 
of its branching ratio at the level of 2 • 1CP 10 would signal the presence of physics beyond 



the Standard Model |52j. The Standard Model sensitivity is expected to be reached at AGS 



around the year 2000 [99|. Also Fermilab with the Main Injector could measure this decay 



The present upper bound on Br^K-^ — » iftvv) is about five orders of magnitude above 
the Standard Model expectation ( 7.37| ). FNAL-E799 expects to reach the accuracy 0(1O~ 8 ) 
and a very interesting new experiment at Brookhaven (BNL E926) [^] expects to reach the 
single event sensitivity 2- 10 -12 allowing a 10% measurement of the expected branching ratio. 
There are furthermore plans to measure this gold-plated decay with comparable sensitivity 



at Fermilab [101] and KEK go 



8 Summary 

We have given a general description of OPE and Renormalization Group techniques as applied 



to weak decays of mesons. Further details can be found in [g]], [2^, 23, |103[] , One of the 
outstanding and important challanges for theorists in this field is a quantitative description 
of non-leptonic meson decays. In the field of K-decays this is in particular the case of 



the AI = 1/2 rule for which some progress has been made in [104]. In the field of B 



decays progress in a quantitative description of two-body decays is very desirable in view of 
forthcoming B-physics experiments at Cornell, SLAC, KEK, DESY, FNAL and later at LHC. 
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Recent reviews on non-leptonic two-body decays are given in [105, 106, 107] where further 
references can be found. 

I would like to thank Peter Breitenlohner, Dieter Maison and Julius Wess for inviting me 
to such a pleasant symposium. 
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Table 1: References to NLO Calculations 



Decay 


Reference 


AF = 1 Decays 


current-current operators 
QCD penguin operators 
electroweak penguin operators 
magnetic penguin operators 
Br(B) SL 

inclusive AS = 1 decays 


H 

© 
a 


ouj, 

1 

41, 


36, 37, H 
36, 37] 

42,|| 


Particle-Antiparticle Mixing 


m 

V3 


m 


47] 


Rare K- and i?-Meson Decays 


K° L -» n°uD, B -» l+r, B -> X B vi> 
K + — * TT^VU, — > 
K + — ► TT+Hp, 

K L -» 7r°e + e- 
B -» 

B^X s7 

inclusive B — > Charmonium 


© 

© 
©- 


1, 

H 

57| 

|65 


51, 52] 


Two-Higgs Doublet Models 


£?^x s7 


© 


6§ 35 1 


Supersymmetry 


AMk and 

s^x s7 


a 



Table 2: Electroweak Two-Loop Calculations 



Decay 


Reference 


k° l -> 7r°w, s -» z+r, B -> JS>z7 

B^X s7 

S° — B° mixing 


|72] 

[73, 74, 75] 
|76] 
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